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HAMMERSLEY'S PROCESS WITH SOURCES AND SINKS 

By Eric Cator and Piet Groeneboom 

Delft University of Technology 

We show that, for a stationary version of Hammersley's process, 
with Poisson "sources" on the positive i-axis, and Poisson "sinks" on 
the positive j/-axis, an isolated second-class particle, located at the 
origin at time zero, moves asymptotically, with probability 1, along 
the characteristic of a conservation equation for Hammersley's pro- 
cess. This allows us to show that Hammersley's process without sinks 
or sources, as defined by Aldous and Diaconis [Probab. Theory Re- 
lated Fields 10 (1995) 199-213] converges locally in distribution to a 
Poisson process, a result first proved in Aldous and Diaconis (1995) 
by using the ergodic decomposition theorem and a construction of 
Hammersley's process as a one-dimensional point process, develop- 
ing as a function of (continuous) time on the whole real line. As a 
corollary we get the result that EL(t,t)/t converges to 2, as t — > oo, 
where L(t,t) is the length of a longest North-East path from (0,0) 
to (t,t). The proofs of these facts need neither the ergodic decompo- 
sition theorem nor the subadditive ergodic theorem. We also prove 
a version of Burke's theorem for the stationary process with sources 
and sinks and briefly discuss the relation of these results with the 
theory of longest increasing subsequences of random permutations. 



1. Introduction. Let L n be the length of a longest increasing subse- 
quence of a random permutation of the numbers l,...,n, for the uniform 
distribution on the set of permutations. As an example, consider the permu- 
tation (5,3,6,2,8,7,1,4,9). Longest increasing subsequences are (3,6,7,9), 
(3,6,8,9), (5,6,7,9) and (5,6,8,9). In this example the length of a longest 
increasing subsequence is equal to 4. 

In Hammersley (1972) a discrete-time interacting particle process was in- 
troduced, which has at the nth step a number of particles equal to the length 
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of a longest increasing subsequence of a (uniform) random permutation of 
length n. This process is defined in the following way. 

Start with zero particles. At each step, let, according to the uniform dis- 
tribution on [0, 1] , a random particle U in [0, 1] appear; simultaneously, let 
the nearest particle (if any) to the right of U disappear. Then, as shown 
in Hammersley (1972), the number of particles after n steps is distributed 
as L n . Hammersley (1972) uses this discrete-time interacting particle pro- 
cess to show that EL n /yfn converges to a finite constant c> 0, which is also 
the limit in probability [and, as noticed later by H. Kesten in his discussion 
of Kingman (1973), the almost sure limit] of L n / y/n. To prove that EL n / yfn 
converges to a finite constant c > is the first part of "Ulam's problem," 
the second part being the determination of c. 

Aldous and Diaconis (1995) introduce a continuous-time version of the 
interacting particle process in Hammersley (1972), letting new particles ap- 
pear according to a Poisson process of rate 1, using the following rule: 

Evolution rule. At times of a Poisson (rate x) process in time, a 
point U is chosen uniformly on [0,x], independent of the past, and the 
particle nearest to the right of U is moved to U, with a new particle created 
at U if no such particle exists in [0,3;]. 

For our purposes the following alternative description is most useful. Start 
with a Poisson point process of intensity 1 on R^_. Now shift the interval [0, x] 
vertically through (a realization of ) this point process, and, each time a point 
is caught, shift to this point the previously caught point that is immediately 
to the right. Let L(x,y) be the number of particles in the interval [0,x] after 
shifting to height y. Then, by Poissonization of the length of the random 
permutation, we get 

L~ =L(x,y), 

where 

— D 

N x ,y = #{points of Poisson point process in [0,x] x [0, y]} = Poisson(xy). 

In an alternative interpretation, L(x,y) is the maximal number of points 
on a North-East path from (0,0) to (x,y) with vertices at the points of the 
Poisson point process in the interior of , where the length of a North-East 
path is defined as the number of vertices it has at the points of the Poisson 
point process in the interior of M% . The reason is that a longest North-East 
path from the origin to (x,y) has to pick up a point from each space-time 
path crossing the rectangle [0,x] x [0,y]. Aldous and Diaconis (1995) call the 
evolving point process y t— > L(-,y), y > 0, of newly caught and shifted points 
Hammersley'' ' s interacting particle process. 
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Fig. 1. Space-time paths of Hammersley's process, contained in [0,x] x [0,y]. 

We can also introduce the evolving point process x i— > L(x, •), x > 0, run- 
ning from left to right. Analogously to the description above of the process 
running up, we shift in this case an interval [0, y] on the y-axis to the right 
through the point process in the interior of the first quadrant, and, each 
time a point is caught, shift to this point the previously caught point that 
is immediately below this point (if there is such a point). By symmetry, it is 
clear that the processes y i— ► L(-,y), y > 0, and x L(x, •), x > 0, have the 
same distribution. 

A picture of the space-time paths corresponding to the permutation 
(5,3,6,2,8,7,1,4,9) is shown in Figure 1. In this case [0,x] x [0,y] con- 
tains nine points, and one can check graphically that there are four longest 
North-East paths (of length 4) from (0,0) to (x,y), corresponding to the 
subsequences (3,6,7,9), (3,6,8,9), (5,6,7,9) and (5,6,8,9). Following a ter- 
minology introduced in Groeneboom (2001), we call the points of the Poisson 
point process in the interior of a -points and the North-East corners of 
the space-time paths of Hammersley's process f3 -points. In fact, the actual 
^-coordinates of the a-points in the picture are different from the numbers 
3, 6, ... , but the ranks of these x-coordinates are given by 3, 6, and so on, if 
we order the a-points according to the second coordinate. 

We use a further extension of Hammersley's interacting particle process, 
where we have not only a Poisson point process in the interior of , but also, 
independently of this Poisson point process, mutually independent Poisson 
point processes on the x- and y-axis. We call the Poisson point process on 
the x-axis a process of "sources," and the Poisson point process on the y- 
axis a process of "sinks." The motivation for this terminology is that we 
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Fig. 2. Space-time paths of Hammersley' s process, with sources and sinks. 



now start the interacting particle process with a nonempty configuration of 
"sources" on the x-axis, which are subjected to the Hammersley's interacting 
particle process in the interior of , and which "escape" through sinks on 
the y-axis, if such a sink appears to the immediate left of a particle (with 
no other particles in between). Figure 2 shows how the space-time paths 
change if we add two sources and three sinks (at particular locations) to the 
configuration in Figure 1. 

The interacting particle process with sources and sinks was studied in 
Section 4 of Groeneboom (2002), where it was proved that, if the intensity 
of the Poisson processes on the x- and y-axes are A and 1/A, respectively, 
and the intensity of the Poisson process in the interior of is 1 , the process 
is stationary in the sense that the crossings of the space-time paths of the 
half- lines M+ x {y} are distributed as a Poisson point process of intensity A, 
for all y > 0. The stationarity of the process was proved by an infinitesimal 
generator argument. It also follows from the computations in the Appendix 
of the present paper. The process is studied from an analytical point of view 
in Baik and Rains (2000) (see Remark 3.1 in Section 3). 

In Section 2 we compare Hammersley's interacting particle process, as 
introduced in Aldous and Diaconis (1995), with the stationary extension of 
this process, with sources on the x-axis, and sinks on the y-axis. However, 
as an intermediate step, we introduce a process with Poisson sources on 
the positive x-axis, but no sinks on the y-axis. From Theorem 2.1 in the 
present paper we can deduce that this particle process, with Poisson sources 
of intensity A on the positive x-axis, but no sinks on the y-axis, behaves 
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Fig. 3. Pat/i of isolated second-class particle in the configuration of Figure 2. 



below an asymptotically linear "wave" of slope A 2 through the /3-points as 
a stationary process. 

In a coupling of the process with the stationary process, having both 
sources and sinks, this wave can be interpreted as the space-time path of 
an isolated second-class (or "ghost") particle with respect to the stationary 
process. For the concept "second-class particle" in the context of totally 
asymmetric simple exclusion processes (TASEP), see, for example, Ferrari 
(1992) or Liggett [(1999), Chapter 3]. The second-class particle jumps to the 
previous position of the particle that exits through the first sink at the time 
of exit, and successively jumps to the previous positions of particles directly 
to the right of it, at times where these particles jump to a position to the left 
of the second-class particle; see Figure 3. The space-time path of the isolated 
second-class particle moves asymptotically, with probability 1, along the 
characteristic of a conservation equation for the stationary process. Here we 
establish a connection with the theory of totally asymmetric simple exclusion 
processes. Although we use similar techniques as used for the study of the 
behavior of second-class particles in TASEP, the situation is in a certain 
sense simpler in our case, since we do not have to condition on having a 
second-class particle at the origin at time zero. 

In a similar way we prove that Hammersley's process, with Poisson sinks 
of intensity 1/A, A > 0, on the positive y-axis, but no sources on the x-axis, 
behaves asymptotically as a stationary process above a wave through the 
/3-points of slope A 2 , if the Poisson sinks on the positive y-axis and the 
points of the Poisson process (of intensity 1) in the interior of are inde- 
pendent. By a coupling argument, these processes can be compared directly 
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to Hammersley's process, as defined in Aldous and Diaconis (1995), which 
has empty configurations on the x- and y-axis. The coupling argument gives 
a direct and "visual" proof of the local convergence of Hammersley's pro- 
cess to a Poisson point process with intensity A, if one moves out along a 
"ray" y = X 2 x, which is the main result Theorem 5 of Aldous and Diaconis 
(1995). The convergence of EL(t,t)/t to 2, as t — > oo, then also easily fol- 
lows. This implies that EL n /y/n converges to 2, a result first proved by 
Logan and Shepp (1977) and Vershik and Kerov (1977). 

In Section 3 we study the /3-points of the stationary Hammersley process. 
For these points we prove a "Burke theorem," showing that these points in- 
herit the Poisson property from the a-points. This allows us to show, using 
a time reversal argument, that in the stationary version of Hammersley's 
process, a longest "weakly" North-East path (allowing horizontal and ver- 
tical pieces along the x- or y-axis) only spends a vanishing fraction of time 
on the x- or y-axis. 

2. Path of an isolated second-class particle and local convergence of Ham- 
mersley's process. Fix A > 0, and let 1 1— > L\(-,t) be Hammersley's process, 
now considered as a one-dimensional point process, developing in time t, 
generated by a Poisson process of sources on the positive x-axis of inten- 
sity A, A > 0, a Poisson process of sinks on the time axis of intensity 1/A 
and a Poisson process of intensity 1 in K+, where the Poisson process on 
the x-axis, the Poisson process on the time axis and the Poisson process 
in the plane are independent. It is helpful to switch from time to time the 
point of view of Hammersley's process as a process of space-time paths in 
]R^_ and Hammersley's process as a one-dimensional point process, devel- 
oping in time. This is somewhat similar to the two ways one can view the 
Brownian sheet. Since the second coordinate can (mostly) be interpreted 
as "time" in the sequel, we will denote this coordinate by t instead of y, 
although, with slight abuse of language, we will continue to call the vertical 
axis the "y-axis," following standard terminology. 

We add an isolated second-class particle to the process, which is located 
at the origin at time zero. A picture of the trajectory of the isolated second- 
class particle for the configuration shown in Figure 2 is shown in Figure 3. 
Theorem 2.1 shows that the space-time path of the second-class particle 
is asymptotically linear with slope A 2 . This is to be expected from results 
on totally asymmetric simple exclusion processes (TASEP), as given in, for 
example, Ferrari (1992). For TASEP Burgers' equation is the relevant conser- 
vation equation in a continuous approximation to the process. The analogue 
of Burgers' equation for a macroscopic approximation to Hammersley's pro- 
cess (with neither sources nor sinks) is 




du(x, t) 



+ u(x,t) 
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where u(x,t) is the intensity of the crossings at (x,t); see Liggett [(1999), 
page 316], where the corresponding equation is given for the integrated in- 
tensity. 

This leads us to expect that, analogously to the TASEP results, 

where Xt is the x-coordinate of the second-class particle, and where de- 
notes almost sure convergence, since in this case the path {(x,t) = (t/X 2 , t) : t > 0} 
is a characteristic for (2.1); compare to, for example, (12.1) in Section 12 of 
Ferrari (1992). 

Theorem 2.1. Let ii— > L\(-,t) be the stationary Hammersley process, 
defined above, with intensities X and 1/A on the x- and y-axis, respectively. 
Let Xt be the x- coordinate of an isolated second-class particle w.r.t. L\ at 
time t, located at the origin at time zero. Then 

(2.2) t^Xt^l/X 2 , t^oo. 

The proof of Theorem 2.1 is based on Lemma 2.1. To formulate this 
lemma we first introduce some notation. Let rjt, t > 0, be the stationary point 
process, obtained by starting with a Poisson point process with intensity 
7 > in (0, oo) at time 0, and letting it develop according to Hammersley's 
process on (0, oo), with Poisson sinks of intensity 1/7 on the y-axis, and 
a Poisson point process of intensity 1 in the interior of the first quadrant. 
Furthermore, let at,t > 0, be the stationary process, coupled to r]t,t > 0, by 
using the same points in the first quadrant as used for rj, and starting with a 
(5/7)- "thickening," 5 > 7, of the Poisson point process with intensity 7 > 
on the x-axis, obtained by adding independently a Poisson point process of 
intensity 5 — 7, and letting at develop according to Hammersley's process 
on (0,oo). To get stationarity for the process a, we replace the sinks on the 
y-axis by a 7/5-thinned set, obtained by keeping each sink with probability 
7/5, independently for each sink. Then the sinks on the y-axis for the process 
a have intensity 1/5. Finally, we let t 1— ► £t be the process of second-class 
particles of rj w.r.t. a, that is, the points of & denote the locations where 
the point process at has extra particles w.r.t. the point process rjt- 

We use the notation r]t[0,x] for the number of particles of rjt in the in- 
terval [0,x] at time t, with the convention that particles, escaping through 
a sink in the time interval [0,t], are located at zero. We define a"t[0,x] simi- 
larly. Furthermore, we use the notation r]t(0,x] (at(0,x]) for the number of 
particles of rjt (at) in the open half-open interval (0, x] at time t. Finally we 
define the "flux" Fe(x,t) of £ through x at time t by 



(2.3) 



F^(x,t) =a t [0,x] -r]t[0,x}. 
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Fig. 4. Processes r\ and £. 



The flux is equal to the number of second-class particles in (0, cc] 

at time t minus the number of removed sinks in the segment {0} x [0,t] 
(through which space-time paths of second-class particles start moving to 
the right). Relation (2.3) is in fact a conservation law. 

A picture of the processes n and £ is shown in Figure 4. In this case the 
process a (inside the rectangle [0, x] x [0,t]) is obtained from the process rj 
by adding two sources at the locations z\ (0) and 22(0) and removing a sink 
at height Sq. The crossings of horizontal lines of the space-time paths of the 
process a are the unions of the crossings of (the same) horizontal lines of 
the space-time paths of the processes 77 and £. 

Lemma 2.1. (i) Let 77 6e Hammer sley'' 's process, defined above, with 
sources of intensity 7 > and sinks of intensity 1/7, and let 8 > 7. VFe 
add independently a Poisson point process of intensity 5 — 7 to the Poisson 
process of sources, and perform a j/5-thinning of the Poisson point process 
of sinks of intensity I/7 on the y-axis. Let a be Hammersley' 1 s process, cou- 
pled to 77, and having the augmented set of sources with intensity 5 and the 
thinned set of sinks with intensity 1/5. Finally, let Z t be, at time t, the loca- 
tion of the second-class particle for which the space-time path starts moving 
to the right through the smallest removed sink. Then 
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(ii) Let rj represent Hammersley' ' s process developing from left to right, 
with sources (on the x-axis) of intensity 7 > and sinks (on the y-axis) of 
intensity 1/7, and let < 5 < 7. We add independently a Poisson point pro- 
cess of intensity 5~ l — 7 -1 to the Poisson process of sinks of intensity 7 -1 , 
and perform a 5 /7 -thinning of the Poisson point process of sources of in- 
tensity 7 on the x-axis. Let a' be the process developing from left to right, 
coupled to rj , and having the augmented set of sinks with intensity S^ 1 as 
sources and the thinned set of sources with intensity 5 as sinks. Finally, let 
Z' t be the location of the second-class particle of a' w.r.t. 7/ ', for which the 
space-time path leaves the x-axis through the smallest removed source (of the 
original process rj). Note that the smallest removed source of r] is a removed 
sink for r( . Then 

lim — = 75 a.s. 

i-^oo t 

Proof, (i) Let x > 0. We have 

r] n [0,nx] 1 
lim = — h cry a.s., 

n — >oo fi 

since r] n [0,nx] equals r] n (0,nx] plus the number of sinks for the process 77, 
contained in {0} x [0, n] (where n is a positive integer), and since rj n (0,nx] 
and the number of sinks contained in {0} x [0,n] have Poisson distributions 
with parameters nxj and re/7, respectively. Here we use the stationarity 
of the process 77, implying that r] n (0,nx] has a Poisson distribution with 
parameter nx-f. Note that, for each e > 0, 
00 

P{\r] n (0, nx] — nxj\ > ne} < 00, 

n=l 

and hence, by the Borel-Cantelli lemma, 

P{\rj n (0,nx] — nx^f \ > ne infinitely often} = 0, 

implying the almost sure convergence of rj n (0,nx]/n to cry, as n — > 00. The 
almost sure convergence to I/7 of the number of sinks for the process 77, 
contained in {0} x [0, n], divided by n, follows in the same way. 
Similarly, 

a n [0,nx] 1 
lim = - + xo a.s. 

n^oo n 

Hence, by (2.3), 

fn ,\ 1. Fc(nx,n) 11 . .( 1 1 

(2.4) lim — ^ = T \-x(5-j) = -(S-j)< — -x } a.s. 

n^oo n 7 I 70 J 
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This limit is negative for < x < 1/(^5) and positive for x > 1/(j5). 

We can number the particles of £ according to their position at time 0, 
so that, for i > 0, particle i is the ith second-class particle to the right of 
the origin at time 0. We then let Zi(t) be the position of the iih second- 
class particle at time t > 0. For i < 0, we let Zi(t), i = 0,-1,-2,..., be 
the second-class particles at time t, for which the space-time paths leave 
the y-axis through the removed sinks So, Si, . . . , respectively, ordering these 
removed sinks according to the height of their location on the y-axis; note 
that Zt = zo(t) (see Figure 4). 

Hence Fe(x,t) has the representation 

(2.5) F^x, t) = #{i>0: Zi {t) < x} - #{i < : Zj(t) > z}. 

Note that second-class particles Zj(-), i < 0, starting their space-time path 
to the right at a removed source in {0} x [0,t], and satisfying Zi(t) G [0, x], 
do not give a contribution to (2.5), since they give a contribution to r]t[0,x] 
as a particle of r)t, located at zero, and a contribution to o"t[0, x] as a particle 
of at in the interval (0,x]. These two contributions cancel in (2.3). It is also 
clear from (2.5) that, for fixed t, the flux F^(x,t) is nondecreasing in x. 

Relation (2.5) shows that F^(Z n ,n) = F^(zo(n),n) is equal to zero at 
each time n, and since F^(nx,n) is nondecreasing in x for fixed n, we get 
from (2.4), 

r Z n 1 

hm — = — - a.s. 

n— >oo n syd 

But, since Zt is nondecreasing in t, we then also have 

(ii) The result is obtained from part (i) by reflecting the processes w.r.t. 
the diagonal, and noting that the reflected processes have the same proba- 
bilistic behavior, but with the role of sources and sinks interchanged. The 
limit l/(7<5) changes to 7$ because of the interchange of x- and y-coordinate. 
□ 

Proof of Theorem 2.1. We couple the process t \-> (L\(-,t),X t ) with 
the process 1 1— > (r)t,crt), where the processes 77 and a are defined as in part (i) 
of Lemma 2.1, and where L\(-,t) = rjt and 5 > 7 = A. Then Zt < Xt, for all 
t > 0, where Z t is defined as in part (i) of Lemma 2.1. This is seen in the 
following way. 

At time zero, we have Zq = Xq = 0. Since the process a is obtained from 
the process 77 by a thinning of the sinks and a "thickening" of the sources, 
and the space-time path of Zt leaves the axis {0} x R + through the smallest 
removed sink, it will leave this axis at a time which is larger than or equal 
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to the time the space-time path of Xt leaves the axis, since the space-time 
path of Xt will leave the axis through the smallest sink in the original set 
of sinks. Note that since a has less sinks and more sources: 

(2.6) ^(0,x]<o- t (0,x], t>0,x>0. 

This means that not only Zt becomes positive at a time that is at least as 
large as the time that Xt becomes positive, but also moves to the right at 
a speed that is not faster than that of Xt. Also note that if Zt jumps to a 
position x > Zt-, an r/-particle jumps over it from a position x' > x. Here 
and in the sequel we use the notation Zt- to denote lim^ Z t ', with a similar 
convention for Xt—. 

If Xt- < x and Zt- < Xt- , Xt will jump to x 1 . Since Zt <x' , Zt can never 
overtake Xt- Note that we can have x' > x if several second-class particles 
are next to each other, without a first-class particle in between. In this case 
Z t does not have to move to the position of the 77 particle, but can move to 
the position of the closest second-class particle to the right of it. 

Hence we have, with probability 1, 

lim inf — - > lim — = — - = — . 

t^oo t t^oo t 70 OA 

Since this is true for any 5 > A, we get 

lim inf — - > — ^ . 

t^oo t X 2 

For the reverse inequality, we switch the role of the sources and the sinks, 
and view Hammersley's process as developing from left to right. This time 
we add independently a Poisson point process of intensity 5' 1 — 7 -1 to the 
Poisson process of sinks of intensity 7 -1 , and perform a 5/7-thinning of the 
Poisson point process of sources of intensity 7 on the x-axis, where 7 = A and 
< 6 < 7, and use the process rf and a', defined in part (ii) of Lemma 2.1. 
Note that rf has the same space-time paths as the process rj, defined above. 
In the coupling we now consider L\ as a process developing from left to right 
and take L\(t, •) = rf t . 

Let X' x be an isolated second-class particle for the process running from 
left to right in the same way as Xt is an isolated second-class particle for the 
process running upward. Trajectories of X and X' are shown in Figure 5. 

We have 

(2.7) X(X'(x))<x, x>0, 

writing temporarily X'{x) instead of X' x and X(u) instead of X u . Equa- 
tion (2.7) is equivalent to noting that the trajectory of (Xt,t) lies above 
the trajectory of (x,X' x ) (see also Figure 5). This follows from the fact that 
if (Xt,t) hits a space-time path at a point North- West of the point where 
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(x,X' x ) hits the same space-time path, this must also be true for the next 
space-time path, since the first trajectory moves up, and the second trajec- 
tory moves to the right. 

By Lemma 2.1 and the argument above, now applied on the process mov- 
ing from left to right, we get the relation 



(2.8 



liminf — > lim — = 6X, 

x — >oo x x — *°° X 



with probability 1. But the almost sure relation liminf^^oo X' x /x > 5X im- 
plies for the process tt—*Xt the almost sure relation 



(2.9) 



X 

lim sup — < 1/(<5A), 

t— >oo t 



since we get for each A' > 1/(5 A), with probability 1, 

X{t/\') X(X'(t)) t 
lim sup j— — < lim sup — — < lim — — 

t^oc t/X' ~ t^oc t/X' -t-wot/A' 



A', 



using (2.8) in the first inequality and (2.7) in the second inequality. 
Since (2.9) is true for any 5 < A, we get, with probability 1, 

lim sup — < — r . 

t^oo t X z 

The result now follows. □ 



(8,1) 
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Fig. 5. Trajectories of (X t ,t) and (x,X' x ). 
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Remark 2.1. The second-class particle X' x , introduced at the end of the 
proof of Theorem 2.1, plays the same role for Hammersley's process, running 
from left to right, as the second-class particle Xt plays for Hammersley's 
process, running up. It therefore has to satisfy 

(2.10) lim ^ = A 2 , 

x— >oo x 

with probability 1. Note that we get an interchange of the x and t coordinate 
which leads to A 2 in (2.10) instead of the 1/A 2 in (2.2), but that the line along 
which (x,X' x ) tends to oo is in fact the same as the line along which (Xt,t) 
tends to oo. 



The following lemma will allow us to show that Theorem 2.1 implies both 
the local convergence of Hammersley's process to a Poisson process and 
the relation c = 2 [which is the central result Theorem 5 on page 204 in 
Aldous and Diaconis (1995)]. 

Lemma 2.2. Let L\ be the stationary Hammersley process, defined in 
Theorem 2.1. Furthermore, let L7 y be the process obtained from L\ by omit- 
ting the sinks on the y-axis, and let L~^ x be the process obtained from L\ 
by omitting the sources on the x-axis. L x y is coupled to L\, by using the 
same point process in the interior of R+, and the same set of sources on 
the x-axis, and L^ x is coupled to L\, by using the same point process in the 
interior o/R+, and the same set of sinks on the y-axis. Then: 

(i) The processes L\ and L^ y have the same space-time paths below the 
space-time path 1 1— ► (X t ,t) of the isolated second-class particle X t for the 
process tt-t L\(-,t). 

(ii) The processes L\ and L~^ x have the same space-time paths above the 
space-time path 1 1— > (t,X' t ) of the isolated second-class particle X' t for the 
process t L\(t, ■), running from left to right. 

Proof. Omit the first sink at location y\ on the y-axis. Then the path 
of L\ leaving through (0, yi) is changed to a path traveling up through the 
/3-point with y-coordinate y\ to the right of (0, yi) until it hits the next path 
of the original process. At this level the path of the changed (by omitting the 
smallest sink) process is going to travel to the left, and the next path will go 
up (instead of to the left) through the closest /3-point to the right. And so 
on. The "wave" through the /3-points that is caused by leaving out the first 
sink is in fact the space-time path of the isolated second-class particle Xt 
(see Figure 3). 

We can now repeat the argument for the situation that arises by leaving 
out the second sink. This will lead to a "wave" through /3-points that is going 
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to travel North of the first wave that was caused by leaving out the first sink. 
This wave is the space-time path of an isolated second-class particle in the 
new situation, where the first sink is removed. Below the first wave the 
space-time paths remain unchanged. The argument runs the same for all 
the remaining sinks. 

(ii) The argument is completely similar, but now applies to the process 
running from left to right instead of up (see the end of the proof of Theo- 
rem 2.1). □ 

In the proof of Corollary 2.1 we will need the concept of a "weakly North- 
East path," a concept also used in Baik and Rains (2000). 

Definition 2.1. In the stationary version of Hammersley's process, a 
weakly North-East path is a North-East path that is allowed to pick up 
points from either the Poisson process on the a>axis or the Poisson process 
on the y-axis before going strictly North-East, picking up points from the 
Poisson point process in the interior R+ . The length of a weakly North-East 
path from (0,0) to (x,t) is the number of points of the Poisson processes 
on the axes and the interior of M. 2 , on this path from (0,0) and (x, t). A 
strictly North-East path is a path that has no vertical or horizontal pieces 
(and hence no points from the axes). 

Note that the length of a longest weakly North-East path from (0, 0) 
to (x,t) in the stationary version of Hammersley's process is equal to the 
number of space-time paths intersecting [0,x] x [0,£], just as in the case of 
Hammersley's process without sources or sinks (in which case only strictly 
North-East paths are possible). 

Corollary 2.1 [Theorem 5 of Aldous and Diaconis (1995)]. Let L be 
Hammer sley" ' s process on R_|_, started from the empty configuration on the 
axes. Then: 

(i) For each fixed a > 0, the random particle configuration with counting 
process 

yi->L(t + y,at)-L(t,at), y>-t, 

converges in distribution, as t — > oo, to a homogeneous Poisson process on R, 
with intensity J~a. 
(ii) 

lim EL(t,t)/t = 2. 

t— >oo 
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Proof, (i) Fix a' > a, and let, for A = \fa' , L7 V be Hammersley's pro- 
cess, starting from Poisson sources of intensity A on the positive rr-axis, and 
running through an independent Poisson process of intensity 1 in the plane 
(without sinks). Then we get from Theorem 2.1 and Lemma 2.2 that the 
counting process y t— > L~^ y (t + y, at) — L^ y (t, at) converges in distribution to 
a Poisson process of intensity A, since the process, restricted to a finite in- 
terval, lies with probability 1 at level t to the right of the space-time path 
of the isolated second-class particle X t , as t — > oo . 

If we couple the original Hammersley process and the process L~^ y via 
the same Poisson point process in the plane, we get that at any level the 
number of crossings of horizontal lines of the process L is contained in the 
set of crossings of these lines of the process L x y , since the latter process has 
sources on the x-axis and no sinks on the y-axis. Hence, for a finite collection 
of disjoint intervals [aj, bi), i = 1, . . . , k, and nonnegative numbers 9\, ... ,9k, 
we obtain 

Eexv\-Y^Oi{L{t + bi,at) - L{t + a h at)}^ 

> £exp j - ]Tj 0, L {LI y (t + bi, at) - L~ y (t + a u at)}| . 
But the right-hand side converges by Theorem 2.1 and Lemma 2.2 to 

^pj-E^-^H 1 - 6 " 91 }}' 

liminf E'exp-! — >^ 6AL(t + bi, at) — L(t + ai,at)\ > 
t— >oo * — ' 

I 4 = 1 ) 



so we get 



(2.11) 



A similar argument, but now comparing the process L with a process L, x , 
having sinks of intensity 1/A = 1 /Va? on the y-axis (which can be considered 
to be "sources" for Hammersley's process, running from left to right), but 
no sources on the x-axis, shows 

limsupi?exp< — Oi{L(t + b{,at) — L(t + Oj, at)} > 

(2.12) 

< e -Eti A ( 6 *-°*){ 1 - e ~ 9i > 
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for any a' < a, since in this case the crossings of horizontal lines of the 
process L are supersets of the crossings of these lines by the process L~^ x . 

That the crossings of horizontal lines of the process L are supersets of the 
crossings of horizontal lines by the process L~^ x can be seen in the following 
way. Proceeding as in the proof of Lemma 2.2, we can, for the process L\, 
omit the sources one by one, starting with the smallest source. The omission 
of the smallest source will generate the path of a second-class particle X[ , and 
the paths of L\ will, at the interior of a vertical segment of the path of X' t , 
have an extra crossing of horizontal lines w.r.t. the paths of the process with 
the omitted source. On the other hand, the process with the omitted source 
will have extra crossings of vertical lines, since some particles will make 
bigger jumps to the left. We can now repeat the argument by omitting the 
second source, which will lead to a further decrease of crossings of horizontal 
lines, and so on. 

Combining (2.11) and (2.12), we find 

lim £exp{ - Y OALit + b h at) - L(t + en, at)}) = e~ Eti^-^v^C 1 -^ 9 ^ 
I U. J 

and the result follows. 

(ii) Since the length of a longest strictly North-East path is always smaller 
than or equal to the length of a longest weakly North-East path, in the 
situation of a stationary process with Poisson sources on the positive x-axis 
and Poisson sinks on the positive y-axis, both with intensity 1, we must 
have, for each t > 0, 

EL(t,t)/t < 2, 

since the expected length of a longest weakly North-East path from (0, 0) 
to (t,t) is 2t for the stationary process. 

The latter fact was proved in Groeneboom (2002), and comes from the 
simple observation that the length of a longest weakly North-East path from 
(0,0) to (t,t) is equal to the total number of paths crossing {0} x [0,t] and 
[0, t] x {t}. Since the number of crossings of {0} x [0,t] has a Poisson(t) 
distribution by construction, and the number of crossings of [0, t] x {t} also 
has a Poisson (t) distribution, this time by the stationarity of the process L\, 
where A = 1 in the present case, we get that the expectation of the total 
number of crossings of the left and upper edge is exactly 2t. 

To prove conversely that limnuv+oo EL(t, t)/t > 2, we first note that 
L(t,t) is in fact the number of crossings of Hammersley's space-time paths 
with the line segment [0,t] x {£}. Take a partition 0,t/k,2t/k, . . . ,t of the 
interval [0,t], for some integer k > 0. Then the crossings of the space-time 
paths of L of the segment [(i — l)t/k,it/k] x {t} contain the crossings of 
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this line segment by the paths of a Hammersley process L7 X with sinks of 
intensity 1/Aj = 1/ y/al, ai < k/i, on the y-axis, but no sources on the x-axis. 

But, by Theorem 2.1 and Lemma 2.2, the crossings of the process L7 X 
with the segment [(i — l)t/k,it/k] x {t} belong, as t — > oo, to the stationary 
part of the process with probability 1, since ai < k/i. 

We now have 

by uniform integrability of t~ 1 L^ x ('yt,t), 7 6 (0,i/fc], t > 0, using, for ex- 
ample, the fact that the second moments are bounded above by the second 
moments of the corresponding stationary process with sources of intensity Aj 
and sinks of intensity 1/Aj. Hence we get, by summing over the intervals of 
the partition, 

1 k 

lim inf EL(t,t)/t > - V Ja~. 
t->oo k 

Letting ai\ k/i, we obtain (still for fixed k) 

k 

lim inf EL(t,t)/t > Vl/V5 = 2(1 + OCL/k)), 

t^oo * — ' 
i=l 

and (ii) follows by letting k — > 00 in the latter relation. □ 

3. Burke's theorem for Hammersley's process. In this section we show 
that, in the stationary version of Hammersley's process with sources on the 
x-axis and sinks on the y-axis, the /3-points inherit the Poisson property 
from the a-points. One could consider this as a version of Burke's theorem 
for Hammersley's process. Burke's theorem [see Burke (1956)] states that the 
output of a stationary M/M/l queue is Poisson. An interesting generaliza- 
tion of Burke's theorem is discussed in O'Connell and Yor (2002). A version 
of Burke's theorem for totally asymmetric simple exclusion processes is given 
in Ferrari [(1992), Theorem 7.1]. Burke's theorem is essentially based on a 
time-reversibility property and for our result on the /3-points this is also the 
case. Our version of Burke's theorem runs as follows. 

Theorem 3.1. Let L\ be a stationary Hammersley process on [0, T\\ x 
[0,T2], generated by a Poisson process of "sources" of intensity A on the 
positive x-axis, a Poisson process of intensity 1/A of "sinks" on the positive 
y-axis and a Poisson process of intensity 1 in M.^_, where the three Poisson 

processes are independent. Let L^ denote the point process of (3-points in 
[0,Ti] x [0,T2], that is, the North-East corners of the space-time paths of 
the process L\, restricted to [0,71] x [0, T2], L A n the entries of the space-time 
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paths on the East side of [0, Ti] x [0, T2] and L^ ut the exits of the space-time 
paths on the North side. Then is a homogeneous Poisson point process 
with intensity 1 in [0,Ti] x [0, T-j], is a homogeneous Poisson process of 
intensity 1/A and T^ ut is a homogeneous Poisson process of intensity A, and 
all three processes are independent. 

Proof. We define a state space E as the possible finite point configu- 
rations on [0,Ti], so E = \J n < Lo E n , where 

E n = {(xi, ...,x n ):0<xi<---<x n <T 1 } (n > 1) 

and Eq = {0}, the empty configuration. We endow each E n with the usual 
topology, which makes E into a locally compact space. We define a Markov 
process (X t )o<t<T 2 011 E such that X t is the point configuration of the 
Hammersley process L on the line [0, Ti] x {t}. In particular we have that 
Xq is distributed according to a Poisson process with intensity A. From 
the definition of the Hammersley process it is not hard to see that the 
generator G of this Markov process is given by 

Gf{x) = f(Ktx) dt + lf(Cx) -(1 + T^j f(x) 

where / G Cq{E), C corresponds to an exit to the left and TZt corresponds 
to an insertion of a new Poisson point at t, so 

C-E^E-Cx= i ( x 2,---,£n), if xeE n (n>2), 
[0, if xe E UEx, 



and for < t < T- 




• j %i— 1) t, Xi-\-\ , . . . , X n ), 

TZ t ■ E — ► E : 1Z t x = { if Xi-\ < t < Xi (x G E n ), 

.,x n ,t), if x n < t (x G E n ). 

Here we use the convention that xo = 0. To prove that G is indeed the 
generator, we fix / G Cq(E) and x G E and consider the transition operators 

P t f(x) = E(f(X t )\X = x) (i>0). 

We will consider the process for a time interval [0, h\ (h | 0) and call Ah the 
number of Poisson points in the strip [0, Ti] x [0,h] and Sh the number of 
sinks in {0} x [0,h]. Then 

Phf(x) = f(x)P(A h = and S h = 0) 

+ ^ f(TZ t x) dt ■ P{A h = 1 and S h = 0) 
J-i Jo 

+ f(Cx)P(A h = and 5 A = 1) + 0(h 2 ) 
= f(x) (l - T\h - l/i) + /(^x) (it + ^/(£x) + 0(/i 2 ). 
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This shows that for every / € Cq(E) and every x € E, 

d_ 
di 

Since X t is clearly a homogeneous Markov process, we get for t £ [0, T?] 

d 



P t f(x) = Gf(x). 

t=o 



< 31 » 



P.f(x) = GP,f(x). 

t 

Now we note that G is a continuous operator on Cq(E), so e tG exists and is 
also a continuous operator. Since 

d 



ds 



e sG f(x) = Ge tG f(x), 

s=t 



(3.1) together with the uniqueness of solutions of a differential equation 
proves that 

P t f(x) = e tG f(x). 

The key idea to prove the theorem is to consider the time-reversed process 
X s = \miX T2 „ s , (X T2 =X ). 

s'ls 

We take the left-limit of the original process X to ensure the cadlag prop- 
erty of (X s )o< s <t 2 - Since, given Xt, the past of the process X is independent 
of the future, it follows immediately that X is a Markov process, possibly 
inhomogeneous. However, if we define /i as the probability measure on E 
induced by a Poisson process of intensity A, then Xq ~ /i and fj, is a station- 
ary measure for the generator G, which implies that X also is stationary 
and homogeneous. The stationarity of X was shown in Groeneboom (2002), 
but will also be a consequence of calculations done in the Appendix. Now 
consider the transition operators 

P t f(x) = E(f(Xt)\Xo=x) (t>0) 

for the time-reversed process. Then, for f,g& Cq(E) and h > 0, 

E(f(X t+h )g(X t )) = E{g(X t )E(f(X t+h )\X t )) 

= E(P h f(X t )g(X t )) 

Phf(x)g(x)/J,{dx). 

E 

We also have 

E(f(X t+h )g(X t )) = E(f(X t+h )E(g(X t )\X t+h )) 
= E(f(X t+h )P h g(X t+h )) 

f(x)P h g(x)n(dx). 
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We use that, due to the stationarity of the process X, Xf and Xt+h both 
have marginal distribution /x. Combining these results gives 

(3.2) f P h f(x)g(x) i x{dx) = [ f(x)P h g(x)fi(dx). 

JE JE 

In the Appendix we calculate the operator G* , defined by the equation 

(3.3) / Gf(x)g(x)fi(dx) = [ f(y)G*g(y)fi(dy) for all f,g£ C (E). 

JE JE 

It is shown there that 

(3.4) G *g(y)= J^ 1 g (C s y)d S + jg(Ky)- (l + T^g(y), 

where in an analogous way as before we define 1Z : E — > E as an exit to the 
right and C s : E — > E as a new point at s such that the point directly to the 
left of s moves to the right. 

We will use (3.4) several times. First of all, since G*l = 0, it shows that 
fx is a stationary measure. Second, we see that for g 6 L°°(/i) 

||G*0||oo <2^ + Ti)||g|| 00 , 

which proves that G is in fact a continuous operator on L 1 ( / u), as well as 
a continuous operator on Cq(E). Since Pt = e tG , Pt is also a continuous 
operator on L 1 ^). Therefore, (3.2) now shows that Pt = P£ = e tG , so in 
fact, using the same argument as before, G = G* . So the reversed process 
has the generator G* . 

Now we define a reflected Hammersley process X v as follows: we take 
the original stationary Hammersley process and reflect all the space-time 
paths with respect to the line segment x [0, T2]; call this a vertical 

reflection. So all points now move to the right and exit on the East side. 
One verifies that the generator for X v is given by G* in the same way we 
did it for the process X, and as X v also starts with a Poisson distribution 
of intensity A, it has the same distribution as X. Note that if one wishes 
to make a picture of the space-time paths of X, one can take the original 
Hammersley process and reflect all the space-time paths with respect to the 
line-segment [0,Ti] x {7^2}, a horizontal reflection. 

Since in X v all the jumps in (0,Ti) x (0,r2) are made toward a point of 
a vertically reflected Poisson process, and in the process X all these jumps 
are made to the horizontally reflected /3-points of the original Hammersley 
process, we have proved that the /3-points are distributed according to a 
Poisson process with intensity 1. Furthermore, in the process X v paths exit 
on the East side according to a Poisson process with intensity 1/A, and this 
corresponds to L A n , horizontally reflected. The process £^ ut , also horizontally 
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reflected, corresponds to the entries of X at the x-axis, and is therefore 
Poisson with intensity A. Finally, the independence of the three processes 
follows from the fact that this is true (by construction) for X v . □ 

Theorem 3.1 allows us to show that a longest weakly North-East path 
from (0,0) to (t/X 2 ,t) only spends a vanishing proportion of time on either 
the x- or y-axis. For the concept of longest weakly North-East path, see 
Definition 2.1. 

Corollary 3.1. Under the same conditions as Theorem 3.1, a longest 
weakly North-East path from (0,0) to (t/X 2 ,t) spends a vanishing proportion 
of time on either the x- or y-axis, in the sense that the maximum distance 
from (0,0) of the point where a longest weakly North-East path leaves the x- 
or y-axis, divided by t, tends to zero with probability 1, as t — > oo. 

Proof. Consider a longest weakly North-East path from (0, 0) to (t/X 2 ,t). 
Such a path can be associated with a path of a second-class particle from 
(t/A 2 ,i) to (0,0) for the time-reversed process, running through the same 
a-points as the longest weakly North-East path, but for which the roles of 
a- and /3-points are interchanged. This means that for the reversed process 
the associated path lies below or coincides with the path of the second- 
class particle that starts moving through the crossing of the upper edge 
[0,t/A 2 ] x {*}, closest to (t/X 2 ,t), moves down to the first a-point on the 
path of the crossing, then moves to the left until it hits the path below the 
highest path crossing the rectangle [0,i/A 2 ] x [0,i], then moves down again, 
and so on. Similarly this path lies above or coincides with the path of the 
second-class particle that starts moving to the left through the crossing of 
the right edge {t/X 2 } x [0,i], closest to (t/X 2 ,t), starts moving down when 
it hits the a-point on the path of the crossing, moves to the left when it hits 
the next path, and so on. 

According to Theorem 2.1 and Remark 2.1, now applied on the reversed 
process, the "/? waves" of the lower and upper path are asymptotically linear 
along the line through the origin with slope A 2 . This implies the statement 
of Corollary 3.1. □ 

Remark 3.1. It is proved inBaik and Rains (2000) that t~ l ^{L\(t, t) - 
2t}, where L\(t,t) is the length of a longest North-East path from (0,0) 
to (t,t) in the stationary Hammersley process (as defined in Theorem 3.1, 
with A = 1), converges in distribution to a distribution function Fq, which is 
related to, but different from the Tracy-Widom distribution function. This 
has the interesting consequence that the correlation between the number of 
points on the left edge and the number of crossings of the upper edge of the 
square [0, t] 2 tends to — 1, as t — > oo. Otherwise the variance of L\(t, t) would 
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be larger than rjt, for some rj > 0, instead of being of order 0(t 2//3 ). We do 
not need their result in our argument, however. Baik and Rains (2000) use 
an analytical approach, applying the Deift-Zhou steepest descent method 
to an appropriate Riemann-Hilbert problem (after using a representation of 
the distribution function in terms of Toeplitz determinants) . This approach 
is rather different from the approach taken here. 

As noted in Baik and Rains (2000), the stationary process is a transition 
between two situations: if the intensities of the Poisson processes on the x- 
axis and y-axis are strictly smaller than 1, we get that i _1 / 3 {L A (i,i) - 2t} 
converges in distribution to the Tracy- Widom distribution. On the other 
hand, if one of these intensities is bigger than 1 (but the intensities are not 
equal), we get convergence of L\(t,t) to a normal distribution, with the 
usual t" 1 / 2 scaling (and a different centering constant). 

Remark 3.2. In Groeneboom (2001) a signed measure process Vt was 
introduced, counting a- and /3-points contained in regions of the plane. The 
Vf-measure of a rectangle [Q,x] x [0,y] is defined as the number of a-points 
minus the number of /3-points in the rectangle [Q,tx] x [0,£y], divided by t. 
The Vj-process has the property that 

V t (S)^V(S), 

almost surely, for rectangles S in the plane, where V is a positive measure 
with density 

(3.5) Mxiy) M^y( Xiy) = _^_, x ' y>0 - 

Here we use the notation V(x,y) to denote the ^/-measure of the rectangle 
[0,a;] x [0, y]. Likewise we write Vt{x,y) for the Vt-measure of the rectangle 
[0,x] x [0,y). 

The problem of proving part (ii) of Corollary 2.1 of the present paper was 
reduced to showing that 



(3.6) / V t (u,v)dV t (u,v) ^ f V(u,v)dV(u,>, 

JB JB 



\c 2 xy, 



where 



V t (u,v)= dV t (u',v'). 

J[0,u]x[0,v) 

Although (3.6) indeed has to hold, the argument for it, given in Groeneboom 
(2001), is incomplete, and needs a result like Theorem 2.1 of the present 
paper to be completed. [The difficulty is caused by the locally unbounded 
variation of the measure Vt, as t — > oo , which has to be treated carefully to 
explain why we need Vt as integrand in the integral in the left-hand side 
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of (3.6) instead of, e.g., Vt, which leads to an integral that is asymptotically 
twice as large.] But since Theorem 2.1 allows us to prove both the local 
convergence to a Poisson process and convergence of EL(t,t)/t to 2, we did 
not pursue the approach in Groeneboom (2001) any further in the present 
paper. 

APPENDIX 

The purpose of this Appendix is to prove (3.4). Remember that 

oo 

E = □ E n 

n=0 

where Eq = {0} and 

E n = {(xi, ...,x n ):0<xi<---<x n < Ti}. 

A Poisson process of intensity A induces a probability measure f/, on E. De- 
note by [i n the restriction of jj, to E n , so jj, n (dx) = X n e~ aTl dx. The generator 
was given by 

G : C (E) - C (E) : Gf(x) = £ f(K t x) dt + ^f(Cx) - Q + Tij f(x). 
Define G + f = Gf + (1/A + Ti)/; we will calculate the dual of G + . Let 

f,geC (E): 

G+f{x)g(x)fi(dx) 

oo „ 

= e- XT ^G + f(0)g(0) + J2 / G + f(x)g(x)» n (dx) 
= e" ATl jf(0)g(0) + e" ATl f(t)g(0) dt 

+ e " ATl £ A ™/ f 1 f{n t x)g{x)dtdx + \ n ~ l f f(Cx)g(x)dx 
= e" ATl \f{0)g{0) + e" ATl f(t)g{0) dt 

oo n ,, 
n= l i= l J{x£E„,x i - 1 <t<x i } 

x g{x) dx dt 

oo „ 

+ e~ ATl A ™ / f(x 1 ,...,x n ,t)g(x)dxdt 

1 J \x£E„.,t>x„\ 



71=1 
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oo „ 

+ e~ ATl £ A"" 1 / f(x 2 , x n )g(x) dx. 

n=l JE " 

Now we make a change of variable for each term in such a way that we 
get f(y) in each of the integrals: 

G+f(x)g(x)fi(dx) 
= e~ Ari jf(0)g(0) + e" ATl £ f(y)g(0) dy 

oo n ,. 

+ e- ATl EE Art / f(y)g(y 1 ,..., m . 1 ,s, 

n= H=l •> {y&E n ,yi<s<y i+1 } 

Vi+i, ■ ■ ■ ,y n ) dy ds 

oo „ 

+ e- ATl EA n / f(y)g(yi,...,y n )dy 

n=l jE n+l 
oo „ 

+ e- ATl E A "" 1 / f(y)g(s, yi ,...,y n . l )dyds 



n=l 



jf(0)g(0)fi o (E o ) + jJ e f(y)g(0)Mdy) 

oo n „ 

+ f(y)g(yi,---,yi-i,s, 

n =i i=i J {y& E n,yi<s<yi+i} 

yi+i,---,yn)Vn(dy) ds 



oo . 

+ E/ f(y)9(s,yi,...,y n )nn(dy)ds 
n=0 J {yeEn,s<yi} 

00 1 f 

+ Et/ f(.y)9(yi,--->yn-i)fj>n(dy) 

oo , Tl s oo , 

= £/ f(y)[ g{C s y)ds)^ n {d y ) + Y^T f{y)g{n y )^ n {dy) 

= J e f(y) g(C*y) ds + jg(Ky)j n(dy). 

Here we define 1Z as an exit to the right and C s as a new point at s such 
that the point directly to the left of s moves to the right, that is, 



1Z:E^E: Tlx 



(xi,. . . ,x n -i), if x £ E n (n > 2) 
0, if x £ EqU E\, 
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and for < s < T\ , 

{{x\ , . . . , Xi—\ , S, Xi-\-\ , . . . , Xn) , 
if Xi < s < Xi + i (x G E n ), 
(s,xi,...,x n ), if s < xx (x G E n ). 

Since G*g = G\g — (1/A + Ti)g, we have shown that 

G*g(y) = g(C s y) ds + ±g(Ky) - Q + T^g(y). 
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